Thermodynamics of Ideal Gas in Cosmology 
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The equation of state and the state functions for the gravitational source are necessary 
conditions for solving cosmological model and stellar structure. The usual treatments are di- 
rectly based on the laws of thermodynamics, and the physical meanings of some concepts are 
obscure. This letter show that, we can actually derive all explicit fundamental state functions 
for the ideal gas in the context of cosmology via rigorous dynamical and statistical calcu- 
lation. These relations have clear physical meanings, and are valid in both non-relativistic 
and ultra-relativistic cases. Some features of the equation of state are important for a stable 
structure of a star with huge mass. 
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The equation of state and the relations between thermodynamic variables are necessary con- 
ditions for solving cosmological model and stellar structure. Obviously to get the reahstic state 
functions for all kinds of matter is a difficult problem due to the lack of the knowledge for some 
materials such as the dark matter. However, for the most important gravitational source in general 
relativity, i.e., the ideal gases or the perfect fluid model, we can rigorously solve the problem. Some 
ordinary but important features of the EOS may be valid for all kinds of matter. 

The cooling mechanism of the expanding universe is one of the most interesting problem of 
the students. The answers given by the standard cosmological text books [l] and the pedagogical 
articles [2, 3, 4] are usually T oc a~^, which is derived from the classical thermodynamics. In 
[5, 6, 7, 8, 9], the author solved the thermodynamical relations according to the Gibbs' law 

TdS = d{pV) + PdV. (1) 

The problem looks underdetermined due to (1) including a number of undetermined quantities. So 
we have to introduce some new concepts and relations such as multi fluid with energy exchange, 
the apparent horizon entropy, the decay of vacuum, to make the equations closed. 

As pointed out in [lU], these approaches seem not convincing the students and satisfying their 
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curiosity. Dividing the velocity of particles into Hubble velocity = ra'{t) and peculiar velocity 
Vpec = <ir'{t), and then analyzing the geodesic of a particle, in [10] the author qualitatively reached 
the conclusion p oc a^^. Then by mass-energy relation E'^ = + m?, he concluded that the cosmic 
temperature should be T oc a^^ for the ultra-relativistic gas but T oc a^'^ for the non-relativistic 
gas. In what follows, we rigorously derive this relation T = T{a) and some other state functions 
of the ideal gases. 

In the Friedman-Robertson- Walker (FRW) space-time, we have the line element 

ds'^ = dt^ - a^{t) {dr^ + fir)de^ + f{r) sin^ Od^f^) , (2) 

where 



f 



sinr if k = 1, 

r if K = 0, (3) 
sinhr if k = —1. 



The evolution equation is given by the Friedman equation 

d'^ + K- ^Aa^ = ^pa^, (4) 
where the over dot stands for ^. The energy- momentum tensor of the perfect fluid and radiation 

Tf"" = {p + P)Uf'U'' - Pgf"" (5) 

satisfies the conservation law 

^ = -3Pa\ (6) 
da 

For a given specific equation of state P = P{p), we can solve the function p = p{a), e.g. for 
radiation P = ^p, then we have p = aa~^. Adding initial values a(to), (4) forms a problem of 
determining the solution. The problem is then thoroughly solved in theoretical sense. 

In the microscopic view, for the ideal gases and photons, the particles are only driven by 
average gravity and move along geodesies, and the collisions between particles can be treated as 
instantaneous behavior. So all thermodynamic functions can be rigorously solved according to 
dynamics and statistics. For convenience, we discuss the thermodynamics of the ideal gases in 
the context of cosmology. If we take (4) and (5) as the macroscopic definitions of p and P, and 
comparing it with the expectation value of the microscopic energy-momentum tensor[l, 11], we 
can inversely derive the equation of state P = P{p) and some other relations. The final results are 
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independent of any concrete space-time, so they are actually valid for the ideal gas in any curved 
space-time due to the principle of equivalence. This is the basic idea of this study. 

At first, we discuss the Gibbs' law (1) in context of cosmology. We have the following conclusion. 

Theorem 1 For the materials with energy-momentum tensor (5), the cosmological principle or 
the FRW metric implies the isentropic process dS = 0. 

This can be checked as follows. Let V = Oa'^, where Q. is any given constant comoving volume 
independent of o, substituting it into (1) and using energy-momentum conservation law (6), we 
find dS = 0. 

To analyze the dynamics of the particles, we introduce a useful lemma 
Lemma 2 If the line element of the orthogonal subspace has the following form, 

ds^ = A{t)dt'^ + g^^{t)dxf'dx'', (7) 

where A and Tj^y only depend on the coordinate t, then the geodesic in this subspace can be solved 
by 



where are constants, and g^'^gua = <J« • 

Lemma 2 can be checked directly. For the FRW metric (2), the line element in the orthogonal 
subspace (t,r) is given by ds'^ = dt^ — a{t)^dr^. According to Lemma 2, we have the geodesic 
equation 

ds ds a 

where C is a constant only depends on the initial data. By (9) we get the drifting speed of a 
particle 

o-dr C ^ vn , , 

v= — = , C = , " g to . (10) 

dt Va^ + C^ 

So the momentum of a particle p = J^"^^2 satisfies 

p{t)a{t)=p{to)a{to), (11) 

where m„ is the proper mass of the particle. For the massless photons, we can check that the 
wavelength X{t) satisfies = so their momentum p also satisfy (11). Although (11) is 
derived in subspacetime {t,r), but it is suitable for all particles due to the symmetry of the FRW 
metric. 
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The relation between momentum p and the kinetic energy K is given by 

p^ = K{K + 2m). (12) 

For fermions, we have the energy distribution of Fermi-Dirac statistics due to the Pauh exclu- 
sion principle. The Fermi-Dirac distribution achieved great successes in the condensed matter 
physics[12], where the influence of Pauli exclusion principle is important. However for the ideal 
gas, the situation is different, where one can hardly connect the drifting movement of the particles 
with Pauli principle. Secondly, the proper mass of a species of particles is a constant, which seems 
inadequate to become a part of random variable, namely the random variable should be kinetic 
energy K rather the total energy E = K + m. Accordingly, the statistical distribution of the kinetic 
energy K for ideal gas should be 



The nonrelativistic limit of (13) is just the Maxwell velocity distribution, in which the random 
variable is the velocity of a particle v. Considering the following reasons, (13) is a more general form. 
First, the partition function of any ensemble is related to the random energy of its components. 
Second, the Maxwell distribution is not Lorentz invariant, but (13) is suitable for both the non- 
relativistic and ultra-relativistic cases. Third, (13) is suitable for multi-particle case. At last, (13) 
is convenient for mathematical treatment. In the following calculation, we only need the first and 
second order moments rather than the concrete J'iK) for the following discussion, so we assume 

K = KT{K)dK = -kT, / K'^r{K)dK = —{kTf, (14) 
Jo 2 Jq 2a 

where cr is a constant depending on specific distribution function !F{K), and for (13) we have 

(7 = |. For the drifting movement of a particle with proper mass m„, by (11) we have Pn = 

where Cn are constants only depending on the initial data at t = io- Then on one hand, for all 

particles we have the mean square momentum directly 

P = (15) 

where Cq is a constant only determined by initial data at Iq. 

One may argue that (15) is probably broken by the collision of the particles. The following 
Lemma shows that (15) holds in statistical sense. 

Lemma 3 The mean square momentum of the ideal gas is independent of the elastic collision 
of the particles. 
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Proof For any elastic collision, we have momentum conservation law 

Pl+P2 = Pl+P2, (16) 

and 

p'i+pl = + Pi + 2(Pi • P2 - Pi ■ P2)- (17) 

Taking average for (17), we have 

p2 = p2 ^ A. (18) 

Since the elastic collision is a reversible process, in statistical sense, we have the exactly equal 
numbers of reversible process of (16), so we also have 

p2 = p2 + A. (19) 

Comparing (18) with (19), we have A = and = P^ . The proof is finished. 
Since collision is finished instantaneously, (15) holds for all time t. 
On the other hand, according to statistical principle, by the moments (14) we have 



/ -J^Kn{Kn + 2mn)HKn)dKn 



where Nn is the number of particles with mass m„ in the volume V = Qa^ , and N = Nn is the 
total number of particles in the volume. Comparing (20) with (15), we get the equation of T{a) as 
follows 

*T(*T + 2.«)=(— ) , (21) 

where m = ^nin is the average mass of all particles, and 6 is a constant only depending on 

initial data. Solving (21), we get 

Theorem 4 The cosmological temperature of universe with ideal gas satisfies 

, ^ arhb'^ , . 

kT = , (22) 

where b is determined by the initial data oq and Tq 



l = J^f2 + ^^ (23) 
ao V crm \ am/ 
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By the theorem we find the cosmic temperature is different from the results simply obtained 
from classical thermodynamics. In what follows we establish the relation between p and T, as well 
as the equation of state. For convenience of the following calculation, we make transformation 
dr = a~^dt to transform the coordinate system into the conformal form, then the line element (2) 
becomes 

= a\T) {dr^ - dr^ - f{r)d9^ - f{r) sin^ Odif'^) . (24) 
The Lagrangian of the cosmological model is given by[l, 13, 14] 

jO- = Y^iR - 2A) -Y.mk^l-vl6\x-Xk), (25) 

k 

where is the coordinate of k-th particle, the scalar curvature 

R = 6^^, (26) 



in which the prime stands for ^ = a^, and the drifting speed of k-th. particle 

/ dr fd9 fsinOdif 



dr ' dr ' dr 



k 



are independent variables related to a for variation. Noticing that ^^{x — X^) oc a ^, by variation 
of / = J Ca^dtdVL with respect to a, we get 

a" + Ka- \ka^ = ^ Yl '"^kJl-vl (27) 

where is any given comoving volume with volume element dil = sin OdrdOdtp, which is inde- 
pendent of a. Substituting (10) into (27), we get 

Multiply (28) by o' and integrate it, again by (10) we have 



= ^ E (^fc + "^fc)a + Ci, (29) 

where C\ is a constant, for classical particles C\ = 0. But for nonlinear spinors Ci < 0, which 
resist the universe to become singular[14]. Making statistical average and using (22), we get 
/2 2 1 A 4 _ 4 _ 8^*^ - 3 0-62 \ ^ 



°" + - jAa. = —pa' = ^1 + 5;;^^^^^ I (30) 
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where p and p are respectively kinematic and static mass density defined by 

P=-^ 5^ ^fc = -^ 5^ (i^fc + mfe), P=-^2^"ifc = ^, (31) 

in which q is the comoving density independent of a. Comparing (30) with (4), we get 
Theorem 5 For the ideal gas in cosmology, the mass density satisfies 

/ 3 ab^ \ / 3kT\ / K 

p = p[l + ^==- ]=p(l + -—]=p(l + - 



ni 



, 3 



3 A;r\ 



= go [kT{kT + 2am)] 2 + - — J , (32) 

where qq is a constant. 

Substituting (32) into (6), we get 

Theorem 6 The equation of state for ideal gas in context of cosmology is given by 

P- -P^' _NkT( kT \ 

2aV^^TF V V 2{am + kT)J' ^^^^ 

^ Qo[kT{kT + 2am)f2 f _ kT \ 

m{kT + 2am) \ 2{am + kT) ) ' ^ ' 

The final expressions of state functions (32) and (34) are independent of any parameters of 
the space-time, and the metric (24) is only used as 'piston-cylinder' to drive the particles in the 
derivation. So the state functions (32) and (34) are actually valid for all cases of ideal gases in 
general relativity. From the above calculation, we have the following conclusions. 

1. In cosmology, the thermodynamic process is isentropic, it is always at the state of maximum 
entropy. 

2. The temperature of the Universe varies as (22), it does not simply decrease in proportion to 
the inverse of scale factor. 

3. The equation of state (34) slightly differs from the usual result of thermodynamics. By (32) 
and (34), we find the EOS P = P{p) satisfying the following increasing and causal conditions 
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\p, (T^oo), 



in which Cs = \J P'ip) is the velocity of sound in the materials, so (32) and (34) are reasonable 
in physics. 



8 



The state functions (32) and (34) can be derived in different ways[]l], and such fluid can form 
stable stars with any given large gravitational mass. (35) is a quite common requirement for the 
EOS of gravitational source. More detailed analysis shows that, for any materials with the EOS 
satisfies the conditions (35), they can always form static and stable star. This conclusion can be 
easily checked as follows. 

For the static spherically symmetrical space-time 

ds^ = B{r)dt 2 - A{r)dr'^ - r^{de^ + sin^ Od^^), (36) 

the exact solutions can be formally solved[15, 16, 17, 18]. In these solutions we usually expressed 
the mass density p and pressure P as the functions of {A,B) and their derivatives. Theoretically, 
such solutions form the general solution to the problem. However, in such expressions the properties 
of the EOS is obscure, and then most of the solutions are unrealistic in physics. 

The realistic static asymptotically flat space-time with spherical symmetry can be converted 
into the following initial problem of an ordinary differential equation system, which can be easily 
solved numerically. If the EOS of the gravitational source satisfies (35), the dynamics for the 
space-time is equivalent to the following initial problem[l, 19], 

M'(r) = AnGpr^, M(0) = 0, (37) 
, (p + P)UTTGPr^ + M) 

= - CI(r-2M)r ' = 

where M{r) is the ADM mass within the ball of radial coordinate r. For any given po > we get 
a unique solution with finite radius R < oo. The metric components are then calculated by 

In the region r > R, we have p{r) = P{r) = 0. Similar to the case calculated in [Ifl], all solutions 
of (37) and (38) are singularity-free. That is, we always have < p < po £ind Rs = 2M{R) < R. 

The condition 7 > 1 is necessary and interesting, which leads to the finite radius of the star 
R < 00 due to Cs{p) 0,{p ^ 0) in (38). But in the case 7 = 1, we get R 00, and then the 
space-time is no longer asymptotically flat. This result reveals that, if EOS satisfies (35), for a star 
with any given large mass, we have regular solutions in equilibrium with finite po- The EOS of the 
matter is decisive for the fate of a star. 
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